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STOCHASTIC SPATIAL MODELS OF HOST-PATHOGEN AND 
HOST-MUTUALIST INTERACTIONS I 

By N. Lanchier and C. Neuhauser^ 
CNRS-Universite de Rouen and University of Minnesota 

Mutualists and pathogens, collectively called symbionts, are ubiq- 
uitous in plant communities. While some symbionts are highly host- 
specific, others associate with multiple hosts. The outcomes of mul- 
tispecies host-symbiont interactions with different degrees of speci- 
ficity are difficult to predict at this point due to a lack of a general 
conceptual framework. Complicating our predictive power is the fact 
that plant populations are spatially explicit, and we know from past 
research that explicit space can profoundly alter plant-plant inter- 
actions. We introduce a spatially explicit, stochastic model to inves- 
tigate the role of explicit space and host-specificity in multispecies 
host-symbiont interactions. We find that in our model, pathogens 
can significantly alter the spatial structure of plant communities, pro- 
moting coexistence, whereas mutualists appear to have only a limited 
effect. Effects are more pronounced the more host-specific symbionts 
are. 

1. Introduction. The diversity and structure of plant communities are 
largely determined by nutrient availability, competition among plants, her- 
bivory and associations between plants and their symbionts. The first three 
have been the focus of much ecological research, both empirical and theo- 
retical. The role of symbionts on diversity and structure has received less 
attention [1] and will be the topic of this paper. 

A symbiont is an organism that lives in close association with a host. It 
can have either a harmful effect, in which case we call it a pathogen, or a ben- 
eficial effect, in which case we call it a mutualist. Symbionts are ubiquitous. 
For instance, more than 90% of terrestrial plants associate with mycorrhizal 
fungi [18], a beneficial association that supplies nutrients to the plant and, in 
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return, carbon to the fungal partner. Fecundity and viability of virtually all 
plants are affected by pathogens, sometimes with devastating effects, such 
as Dutch elm disease or chestnut blight that point to the important role of 
pathogens in structuring plant communities. Mathematical models play an 
important role in elucidating the roles of symbionts in community dynamics. 

Modeling of disease dynamics has had a long tradition, starting with the 
model by Kermack and McKendrick [14], which describes the course of a dis- 
ease outbreak caused by a single disease infecting a single host. This model 
and its extensions have yielded enormously valuable insights into disease 
dynamics and potential control strategies. Although originally developed to 
describe epidemics in human populations, it can equally well be applied to 
plant diseases. A key concept of disease dynamics is the basic reproductive 
rate Rq, which is defined as the expected number of secondary infections 
caused by an infected individual when introduced into a population of sus- 
ceptible individuals [2]. The condition for a disease outbreak is given in 
the biological literature as i?o > 1- This condition is based on a single-host, 
single-disease model in a nonspatial population. Including spatial structure 
in the form of local interactions has shown that for a disease to spread, 
Rq needs to exceed a threshold that is greater than that for a nonspatial 
population. The reason for this is the lack of sufficient numbers of suscep- 
tible individuals near the location of a disease outbreak once the disease 
starts spreading. One of the first models where this has been demonstrated 
mathematically is the contact process [11, 16]. 

Much of the theoretical work in the epidemiological literature has focused 
on single-host, single-disease dynamics. A rapidly increasing empirical body 
of work on multispecies host-disease dynamics necessitates the development 
of a theoretical framework. This has only begun recently (see [13], and ref- 
erences therein). An attempt for a broad classification was made by Holt et 
al. [13] using a graphical isocline framework that allows for generalizations 
that are applicable to a wide range of host-pathogen models. A recent model 
by Dobson [4] investigates persistence of a pathogen that can infect multiple 
hosts. 

Modeling of host-mutualist interactions has primarily focused on evolu- 
tionary questions, such as the evolution of cheaters (i.e., symbionts that 
receive benefits but do not confer them). Almost no modeling has been 
done on the effects of local (spatial) interactions on host-mutualist dynam- 
ics. Similarly, no theoretical framework has been developed for multispecies 
host-mutualist interactions. 

An important component of multihost, multisymbiont models is the de- 
gree to which different symbionts and hosts can associate with each other. 
This is referred to as specificity. A specialist symbiont associates with a very 
small number of hosts; a generalist symbiont associates with many hosts. 
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The ease of transmission of a symbiont to a host, referred to as transmissi- 
bility, is another important factor in host-symbiont interactions. 

In this paper we investigate the role of spatial structure caused by local in- 
teractions, such as symbiont transmission and host dispersal, on persistence 
of host-symbiont associations for both generalists and specialists in mul- 
tihost, multisymbiont models. We employ the simplest of all multispecies 
models to describe the host dynamics, the voter model [3, 12]. The voter 
model is defined on the d-dimensional integer lattice, where each lattice 
site is occupied by an individual characterized by one of a finite number 
of types. Individuals give birth to offspring of their own kind at a constant 
rate, and their offspring displace randomly chosen individuals within their 
dispersal neighborhood. The dynamics imply that all sites remain occupied 
at all times. Into this population, we introduce symbionts with varying de- 
grees of specificity and transmissibility. More precisely, our spatial model 
is a continuous-time Markov process '■ — > {1, 2, . . . , k} x {0, 1, ... , k} 
where the integer k denotes both the number of hosts and the number of 
symbionts involved in the interaction. A site x € Z'^ is said to be occupied by 
an unassociated host of type i, i = 1, 2, . . . , k, if ^{x) = (i, 0), and by a host 
of type i, i = 1, 2, . . . , K, associated with a symbiont of type j, j = 1, 2, . . . , k, 
if ^(x) = (i, j). We will write ^f(x) = {(^j (x) , S,'^ (x)) , where Cl{x) denotes the 
type of the host present at x at time t and S,t{x) denotes the type of the 
symbiont present at x at time t, with Cf{x) = denoting the absence of 
a symbiont. We set ||x|| = supj^j^ 2,...,d l^il- The evolution at site x S Z'^ is 
described by the transition rates 

(i,j)^(/c,0) at rate A ^ ^ l{«W=(fc,o)} + ^{«W=(fc/)} f 

0<\\x-z\\<Ri I £=1 ) 

K 

{i,0)^{i,j) at rate aj ^ H ^{?(^)=(Aj)}- 

0<\\x-z\\<R2 £=l 

The transition {k,0) is the birth of an unassociated host at x by 

either unassociated or associated neighboring hosts. The birth rate of unas- 
sociated hosts is equal to A. The parameter g indicates the variation of the 
birth rate of hosts associated with a symbiont. If < 5 < 1, the symbiont is 
a pathogen; if g = 1, the symbiont has no effect on the birth rate of the host 
and we refer to this as the neutral case; if 5 > 1, the symbiont is a mutualist. 
The transition {i, 0) — > is the transmission of a neighboring symbiont j 
to an unassociated host of type i at x. The parameters Cij denote the rate 
at which symbiont j infects host i. This parameter will allow us to mimic 
specialist and generalist symbionts. Births and infections occur within a lo- 
cal neighborhood, with Ri denoting the birth range of hosts, and R2 the 
infection range of symbionts. Neighborhoods are punctured boxes with side 
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2Ri + 1,^ = 1,2, centered at site x, that is, A/"* = {z G Z*^ : < ||x — z\\ < Ri}. 
The cardinahty of this set is denoted by vr^ = 

Before we describe the behavior of the spatially explicit, stochastic model, 
we will look at the mean-field model [9]. The mean-field model is described 
by a system of differential equations for the densities of unassociated and 
associated hosts. To define it, we let Ui denote the density of unassociated 
hosts of type i, i = 1,2, . . . , k, and Vij denote the density of host i associ- 
ated with symbiont j, i,j = l,2,...,K. It follows from the dynamics of the 
spatially explicit, stochastic model that at all times 

i id 

Furthermore, we assume that for i = 1, 2, . . . , k, = /3, and ior i,j = 1,2, n 
with i 7^ j, Cij = a with < a < /?. We define 



u.=^Ui, v.j=^Vij, Vi.=^ Vij and v..=^^Vij. 

i=l i=l j=l i=l j=l 

One way to obtain the mean-field limit is to set the neighborhood ranges 
Ri and R2 equal to R and then let R go to infinity. To obtain a meaningful 
limit, we also need to rescale the parameters A, a and /? by the neighborhood 
size ur; that is, we set = 7^ (this also sets the time scale), and define 

a = — and (3 = — . 

In the limit R ^ 00, sites become independent. If, in addition, the initial 
configuration is translation invariant, the dynamics of the densities for j 
is then described by the following system of differential equations, called 
mean- field equations: 

^ = (1 - Ui){ui+gvi.) - Ui^{uj + gvj.) - bmv.i -a^mv.j. 



dt 

dvii 

— = buiV.i - Vii[u. + gv..), 
dvij 

= auiv.j - Vij{u. +gv..). 

When a = 0, the symbionts are specialists. As a increases to h, the associ- 
ation turns into a generalist relationship. The following results are proved 
in Section 2. When g = l, the system has a conserved quantity, namely the 
initial host densities hi = Ui + Vi., i = 1,2, . . . , k. If (k — l)a + b> k, then for 
g ^ 1, there exists a nontrivial equilibrium with ui = U2 = ■ ■ ■ = > and 
vi. = V2. = ■ ■ ■ = ffi;. > such that for i = 1,2, . . . , n, 

Ui = J and hi = —. 

[K — l)a + b — k{1 — g) k 
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Fig. 1. Solution curves of the mean-field model with a = and & = 3. The parameter g 
is equal to 0.5 and 2, respectively. 

Furthermore, for j 

Vij _ a 
Vii b' 

Numerical simulations indicate that the nontrivial equilibrium is locally sta- 
ble for a <b when 5 < 1, but not for > 1. In addition, if = and {n— l)a-\- 
b> K, all hosts will be associated at equilibrium. If initially only two hosts 
and one symbiont, say symbiont 1, are present, then for < g < 1, simu- 
lations indicate that pathogen 1 will go extinct and both host 1 and host 
2 may coexist. We will conjecture similar behavior for the spatial model, 
except in the one-dimensional, nearest-neighbor case (see Theorem 3 and 
discussion following the theorem). When g > 1, host 2 goes extinct provided 
the initial density of mutualists of type 1 is sufficiently large. The spatial 
analogue of this result is proved in Theorem 2. Both results are illustrated 
in Figure 1. 

We now return to the spatially explicit model with parameters A, a and 
/3. To define the time scale, we set A = 1. We will discuss both the generalist 
case a = (5 and the specialist case a = 0, and provide comparisons with the 
mean-field model. 

Generalist interactions. We consider the generalist case a = /3 of the spa- 
tially explicit, stochastic model. In the corresponding mean- field model, co- 
existence of hosts and symbionts is possible when 6 > 1. For 6 < 1, symbionts 
are unable to persist. When a = b, symbionts are no longer host-specific, and 
the mean-field model can be reduced to one with Ui and Vi.. It is not hard 



6 



N. LANCHIER AND C. NEUHAUSER 



to see then that the nontrivial equihbrium of unassociated hosts, {u.,v..), is 
locally stable since in this case 

^=v..[{b-l)u.-gv..]. 

With u. + v.. = 1, it follows that the boundary point v.. = is unstable for 
b> 1 and that 

^"•'"■^=(6^' 5^) 

is locally stable. Furthermore, any vector (/ii, /i2, . . . , /ir) with hi > and 
J2i hi = l gives rise to an equilibrium if we set 

b-l+g g 

The behavior of the spatially explicit model is more complicated and may 
depend on the spatial dimension, as we will see in the following. 

If g = 1, then the symbionts have no effect on the hosts, which means 
that the spatially explicit processes and are stochastically indepen- 
dent. Moreover, by looking at the transition rates, it is easy to see that 
is a multitype voter model run at rate 1, and that is a multitype con- 
tact process in which particles give birth at rate /Ji^/jj die at rate z//fj . 
See, respectively, [12, 17] for a study of these two processes. It follows that 
there exists a critical value (3c{Ri,R2) G (0, oo) that depends on the neigh- 
borhood sizes z^/jj and z^/jj such that the symbionts can survive if and only 
if /3 > f3c{Ri, R2)- If we ignore host and symbiont types but rather focus on 
associated versus unassociated hosts, then for /3 > Pc{Ri, R2), regardless of 
the spatial dimension, there exists a nontrivial stationary distribution of as- 
sociated and unassociated hosts. Moreover, if d > 3, there exists a stationary 
distribution in which all hosts and symbionts coexist. 

Unfortunately, we cannot say much about coexistence when g ^ 1- To 
analyze this case, we define the "color-blind" process where a site is in state 
if it is occupied by an unassociated host, and in state 1 if it is occupied 
by an associated host. We obtain a particle system Ct'-'^'^ — > {0)1} with 
transitions at x G Z*^: 

— > 1 at rate /3 ^ l^i^(^z)=i}j 

0<\\x-z\\<R2 

1^0 at rate ^ {]l|^(^)=o| c;]l{(;(^)=i}}. 

0<\\x~z\\<Ri 

When g = 0, the process reduces to a biased voter model. When g = 1, it 
reduces to a contact process with birth rate Pi'R2 and death rate i^/j^. We 
denote the critical value of this contact process by (3c{Ri,R2) as above. A 
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standard coupling argument allows us to compare the processes with 5 7^ 1 
and g = I, and to deduce that if g < 1 and /3 > l3c{Ri, R2), then Q has a 
nontrivial stationary measure, while if 5 > 1 and /3 < (3c{Ri,R2), then the 
mutualists die out, that is, Ct^ Sq, the "all 0" configuration [i.e., ("(x) = 0]. 
To cover the remaining cases, we introduce the contact process rjt in which 
particles give birth at rate f3i'R2 and die at rate gi^R^ ■ Then 7]t has a nontrivial 
stationary measure if and only if /3 > gl3c{Ri, R2) which, with a new coupling 
argument, implies that ii g <1 and (3 < g(3c{Ri, R2), then the pathogens die 
out, while if (7 > 1 and /? > g(3c{Ri, R2), then (^t has a nontrivial stationary 
measure. 

We now focus on the case (3 > 1 and 5 > close to 0. First of all, we 
observe that if g = and Ri = R2, then the process Q is the biased voter 
model with parameters /? and 1. In particular, = 0) — > 1 if /3 < 1 while 

P{Q{x) = 1) — > 1 if /? > 1 provided we start with infinitely many O's and I's 
at time 0. Moreover, in the latter case, fixation occurs for the process since 
hosts associated with pathogens are now sterile. The behavior is identical 
to that of the mean-field model. We will use a perturbation argument in 
Section 4 to show that if /3 > 1 and g > is sufficiently close to 0, then 
the pathogens still survive. The results are summarized in Figure 2 and in 
the following theorem where denotes weak convergence and 5q is the 
distribution that concentrates on the "all 0" configuration. 

Theorem 1. Assume that a = f3 and that Co is translation invariant 
with P(Co(2;) = 1) > 0. 

(a) If g <1, then Ct ^ ^0 if P ^ gPc md a nontrivial equilibrium exists 
if (3 > (3c- If g> 1 , then if [3 < [3c, and a nontrivial equilibrium exists 
if (3>g(3c. 

(b) If P > 1, there exists gc > such that if g < gc, then (t ^ with 

m(C(x) = i)/o. 

Part (b) of this theorem will be proved in Section 4. 

Specialist interactions. In the specialist case a = and (3 >0, the process 
is more difficult to investigate since the evolution of each symbiont strongly 
depends on the configuration of the host population. That is, there is no 
particle system Q '■'^'^ — > {0, 1} which allows us to describe the global evo- 
lution of the symbionts regardless of their type. Since for any i = 1, 2, . . . , k 
the symbiont i can live only through hosts of type i, it is, however, easy 
to deduce from a coupling argument that if g <1, then the processes with 
a = and a = (3 can be defined on the same space so that, starting from the 
same configuration, the process with a = has fewer pathogens. In words, 
the survival of the pathogens is harder to obtain with specialist interactions. 
In particular, if < 1 and /3 < max(/3c5, 1), then the pathogens die out. 
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Fig. 2. Phase diagram of the spatial model with generalist interactions. 



The next step is to extend the results of the mean-field model summarized 
in Figure 1 to the corresponding spatial model with short-range interactions. 
To do this, we consider a population of two hosts with only one type of 
symbiont, say symbiont of type 1, and start the evolution with all the hosts 
of type 1 associated with a symbiont. Then, in the limiting case /3 = oo, 
the transition (1,0) — > (1,1) is instantaneous; that is, unassociated hosts of 
type 1 become instantaneously associated with a symbiont, provided that 
Ri < R2 to avoid the problem of an isolated unassociated host that cannot 
be reached by any symbiont. This implies that the process is a biased voter 
model in which hosts of type 1 give birth at rate gi'Rj and hosts of type 2 at 
rate vr^. In particular, if g < 1, then =^ ^2, the "all 2" configuration, while 
if 5 > 1, then ^ 61, the "all 1" configuration. Theorem 2 tells us that, in 
any dimension, the result still holds if g > 1 and /3 E (0, 00) is sufficiently 
large. 

Theorem 2. Assume that a = 0, g > 1 and Ri< R2- At time 0, i^ix) = 
(1,1) or (2,0) for all x G . Then, there is a critical value PjJ^'^ & (0,oo) 
such that if P> PjJ^^, then 61 . 
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The arguments in our proof, however, do not imply the analogous result 
for g < 1. We think actually that except in the one-dimensional nearest- 
neighbor case, = 2) 7^ 1. Before explaining our intuition, we de- 
scribe the behavior of the one-dimensional process with nearest-neighbor 
interactions for the hosts and short-range interactions for the pathogens. 
Assume that at time 0, ^o(^) ~ ^ ^'^^ x > with infinitely many pathogens 
of type 1, and ^q(x) = 2 for x <0 with no associated pathogens. Denote by 
r| = sup{x G = 2} the rightmost host of type 2. Then Tq = 0. The 
following result implies that for /? large enough, 2's will spread to the right 
and eliminate I's together with their associated pathogens. 

Theorem 3. Assume d=l and i2i = i?2 = 1. If a = and g <1, there 
exists Pj^^ S (0, oo) such that if P > Pcr^^ , then — > oo as t —> oo with 
probability 1. 

We conjecture that this result should only be true in the one-dimensional 
nearest-neighbor case. Here is the intuition. We first observe that except in 
the one-dimensional nearest-neighbor case, the dynamics produce isolated 
hosts, that is, hosts of type 1 (resp. 2) surrounded by a cluster of hosts of 
type 2 (resp. 1). As suggested by Theorem 2, when g > I, isolated 2's are 
swallowed very quickly by surrounding I's. On the other hand, when g <1, 
an isolated 1 located in a linearly growing cluster of 2's cannot be invaded 
anymore by any pathogen as soon as the cluster has reached some critical 
size. In this expanding region, the process then behaves like an unbiased 
voter model in which I's can now compete with 2's, and survive. See Figure 
3 in Section 5 for simulations in d = 2. 

Finally, we investigate the coexistence of symbionts in the neutral case 
g = 1- We observe that, in this case, the first coordinate process S^} performs 
a voter model run at rate 1. In particular, in d < 2, there does not exist any 
stationary distribution ^ such that ^{S}{x) = z) 7^ for any z G {1, 2, . . . , k} 
(see [12]). Since in the specialist case, it{x) = i implies that il{x) = i, the 
same conclusion holds for the symbionts. In d>3, coexistence occurs for 
the process , that is, there is a stationary measure fi which satisfies the 
condition above. However, due to the formation of clusters, the problem 
of coexistence of the symbionts remains a difficult question. Namely, the 
voter model exhibits clusters whose diameter can exceed some critical 
size, which prevents the symbionts from spreading out. To get around this 
difficulty, we introduce a modification of the particle system, denoted by 
in which the symbionts evolve as previously but where the hosts now 
perform a threshold 9 voter model. More precisely, the process evolves 
according to the following transitions at a; G Z*^: 




at rate 




if card{2; €Z'^:0<\\x- z\\ < Ri 

and S}{z) = k}>e 



otherwise. 
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(i,0)^('i,i) at rate (3 card{z e Z'^ : < \\x - z\\ < R2 and (z) = i} . 

The introduction of this particle system is motivated by Theorem 1 of [7] 
which imphes that the threshold 9 voter model has a nontrivial stationary 
distribution which is close enough to a product measure to produce our next 
result. We will prove the following result. 

Theorem 4. Let g = 1 and 9 < vr^/ k. If Ri and R2 are sufficiently 

large, then there is a critical value p]!^"^ G (0, 00), depending on k, such that 

the following holds: If /3 > (3^^'^ , then coexistence occurs, and if (3 < I3cr^^ , 
then coexistence is not possible. 

Here coexistence means that there is a stationary measure such that 
fJ-{q{x) = i) ^ for any type i £ {1,2, Unfortunately, we do not 

know how to prove something better than coexistence is not possible when 
/3 < f3j^^- However, we conjecture that starting from product measure in 
which each host is associated with a symbiont and has density 1/k, all the 
symbionts die out. To justify our conjecture, we observe that the processes 
r]l defined by r]l{x) = 1 if S,f{x) = i and r]\{x) = otherwise, do not interact 
since they are confined to their associated hosts. Since Theorem 1 of [7] 
tells us that the hosts coexist with density of each type close to each 
symbiont should remain subcritical. Finally, since the symbionts can only 
spread out through their host, we conjecture that Pj^^ is increasing with 
respect to k. 



Comparison of the spatially explicit and the mean- field model. Numerical 
simulations of the mean-field model indicate that coexistence is only pos- 
sible when g <1. Simulations of the spatially explicit model show similar 
behavior. When g <1 and a < (3, then coexistence of hosts and pathogens is 
possible. We observed that in this case, cluster size is limited by the presence 
of pathogens: In the absence of pathogens, clusters grow at the expense of 
neighboring clusters that contain symbionts. Upon invasion by the preferred 
symbionts (those with infection rate (3), the clusters appear to shrink again. 
The case 5 > 1 and a < (3 \s quite different. Clusters of hosts with their pre- 
ferred mutualists form and appear to continue to grow, just as in the voter 
model case. Less preferred mutualists (those with infection rate a) do not 
seem to be able to persist with preferred mutualists, just as is the mean-field 
case of one host and two symbionts with infection rates a and 6, respectively 
{a <b). In summary, pathogens have the ability to alter the spatial struc- 
ture of their hosts by promoting local diversity, whereas mutualists do not 
alter the spatial structure of their hosts. This difference in behavior is more 
pronounced the more host-specific the symbionts are. 
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The rest of this paper is devoted to proofs. In Section 2 we wih investigate 
the mean- field model. In Section 3 we will prove a preliminary result about 
the biased voter model to prepare the proofs of Theorems 1 and 2 which will 
be carried out in Sections 4 and 5, respectively. Section 6 will be devoted 
to the proof of Theorem 3. Finally, we will investigate the coexistence of 
symbionts and prove Theorem 4 in Section 7. 

2. The mean-field model. The mean-field model was introduced in Sec- 
tion 1. Our first claim was that the host density hi = Ui + vi. is a conserved 
quantity when g = \. K straightforward calculation shows that if g = 1, 

d 

— (Uj +-L'i.) =0, 

from which our claim follows. 

We summarize the behavior of the mean-field model in the following 
proposition. 

Proposition 2.1. For {k — l)a + b> k and g >0, there exists a non- 
trivial equilibrium with ui = U2 = ■ ■ ■ = > and Vi. = V2. = ■ ■ ■ = Vk. > 
such that for i = 1,2, . . . , k, 

Ui = r ; ana hi = —. 

[K — l)a + b — K[l — g) K 



Furthermore, for i^ j 



Vjj _ a 
Vii b ' 



Proof. If we denote by hi = Ui + Vi. the density of host i (both associ- 
ated and unassociated) , then 

dhi 

— =Ui + gvi. - hi{u. + gv..). 
dt 

By setting the right-hand side equal to 0, we obtain 

Ui + gvi. 

hi = . 

u. + gv.. 

It follows that 

hi Ui{l-g)+ghi 



hj Uj{l-g)+ghj 
from which we conclude that 

hi _Ui _ Vi. 
hj Uj Vj. ■ 
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In the symmetric case, hi = h2 = • • • = hf^ = we find ui = U2 = ■ ■ ■ = 
and vi. = V2- = ■ ■ ■ = v^.. The nontrivial equihbrium can then be computed 
expHcitly. We find 

Ui = — and hi = -. 

[K — l)a + — n[l — g) K 

Specifically, when (7 = 0, Uj = and consequently all hosts will be associated 
at equilibrium. The condition for the existence of a nontrivial point equilib- 
rium, namely [n — l)a + b> k, follows directly from requiring that Ui < 1/k 
and hi = Furthermore, it follows from 

buiV.i = Vii{u. + gv..) and auiV.j = Vij{u. + gv..) 

that 

a _ v.i Vij 
b v.j Vii' 

Since v.i = v.j by symmetry, the last claim follows as well. □ 

3. Preliminary results about the biased voter model. As explained in 
the Introduction, if a = /? and Ri = R2, then the "color-blind" process Q 
performs a biased voter model when g = 0. If we set a = and consider a 
population of two host types with only one symbiont type, then the process 

, which describes the evolution of both host types, performs a biased voter 
model in the limiting case (3 = 00. So, to prove Theorems 1 and 2, we will 
start by proving a general result about the biased voter model, and then 
apply a perturbation argument to extend this result to the region g > 
small in the first case, and to the region /3 < 00 large in the second case. 
Let f5i, P2 ^ (0,00), and r]t:'^'^ — > be the biased voter model with 

parameters f3i and (32, that is, the process whose state at site x changes as 
follows: 

i^j at rate /Sj ^ 

0<\\x-z\\<Ri 

It is a well-known fact that if /3i > /32, then P{r]t{x) = 1) ^ 1 as t ^ 00, 
provided that at time 0, the process has infinitely many I's and 2's (see, 
e.g., [6], Chapter 3). 

To prove Theorems 1 and 2, we will follow the strategy described in 
[10], Section 3. We begin with a rescaling argument to estimate the rate 
of convergence of P{r]t{x) = 1). This estimate will have to be good enough 
so that a perturbation argument can be applied. The basic idea is to show 
that for given e > 0, members of the family of processes under consideration, 
when viewed on suitable length and time scales, dominate an M-dependent 
oriented percolation process in which sites are open with probability 1 — 
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e ([8], Section 4). To compare the process with a percolation process, we 
consider a positive integer L to be fixed later, and scale space by setting 

B=[-L,L]'^, ^z)=Lz, B, = '^{z)+B 

for any z G Z"^. Let F be a positive integer, and say that {z,n) is occupied 
if all sites in Bz are occupied by I's at time nTL. The first step in proving 
Theorems 1 and 2 is the following. 

Proposition 3.1. Let e > and Pi > (32- Then M , L and V can he 
chosen in such a way that the set of occupied sites dominates the set of open 
sites in an M -dependent oriented site percolation process where sites are 
open with probability p= 1 — 2e/3. 

The key to the proof is duality ( [6] , Chapter 3) . To define the dual process 
of the biased voter model, we consider two collections of independent Poisson 
processes {T^'^:n > 1} and {U^'^:n > 1} where < ||x — z\\ < Ri, with 
parameter P2 and Pi — /?2, respectively. At times T^'^ we draw an arrow 
from z to X and put a 5 at site x, while at times U^'^ we draw an arrow 
from z to X without putting a 5 at the tip. The process is then obtained 
from the graphical representation as follows: At time T^'^, the particle at 
X imitates the one at z. At time C/^'^, the site x becomes occupied by a 
particle of type 1 if 2; is. We say that there is a path from (x,0) to {z,t) if 
there is a sequence of times sq = < si < • • • < = t and spatial locations 
xq = x,xi, . . . ,Xn = z such that the following two conditions hold: 

1. For i = 1,2, . . . ,n, there is an arrow from to Xj at time Sj. 

2. For i = 0, 1, . . . , n, the vertical segments {xi} x (sj, Sj+i) do not contain 
any 5's. 

Finally, we say that there exists a dual path from (x, t) to {z,t — s), < s <t, 
if there is a path from (z, t — s) to (x, t), and define the dual process starting 
at (x, t) by setting 

= G Z'^ : there is a dual path from (x, t) to {z, t — s)} 

for any < s <t. The reason why we introduce the dual process is that it 
allows us to deduce the state of site x at time t from the configuration at 
earlier times. More precisely, 

rjt(x) = 1 if and only if r]t-s(z) = 1 for some z E Os^'*^. 

See [6], Chapter 3. The strategy to proving Proposition 3.1 can then be 
summarized as follows: Let T = TL and x € B^ with ||2:|| = 1. Then, we will 
prove that, with probability arbitrarily close to 1, there exists a dual path 
As starting at (x,T) and landing in the target set B. More precisely, we 
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will prove that As hits the set J = [— i?i,-Ri] by time T where Ri < L/2, 
and then stays inside B until time T. Recall that i?i denotes the range 
of the interactions. In particular, if B is void of 2's at time 0, then, with 
probability close to 1, B^ will be void of 2's as well T units of time later. To 
define the dual path Ag, we start the process at Aq = (x, T) and go down the 
graphical representation. If Ag comes across a 5 at some time s = T — T^'^ 
with X = As, then move As to z. If As meets the tip of an arrow that is 
without a (5 at some time s = T — U^'^ , then move As to z only if it takes it 
closer to 0. Intuitively, this should cause ^4^ to drift toward the set B. We 
now make this argument precise in a series of lemmas. 

Lemma 3.2. Assume that x £ B^, \\z\\ = 1 and (5i > (32- There exist Ci, 
7i G (0, oo) such that 

sup Px{As ^ J for all s <T) <Ci exp(— 71L) 

for L and V sufficiently large. Here, the subscript x indicates the starting 
point. 

Proof. Let ak denote the A:th time As encounters the tip of an arrow 
(with or without a. 5). At time Cfe, the arrow does not have a 5 at its tip with 
probability {(3i — (32) /Pi > 0. Moreover, if A^i^ ^ J and the arrow does not 
have a d at its tip, then with probability at least l/2d > 0, ^4^ moves closer 
to 0. In particular, if = infj/c > 1 : ylg-fc £ J}, then there is c > such that 

P(iV>cL) <C2exp(-72L) 

for suitable 6*2,72 G (0,oo). Since P{ak — c^k-i > t) = exp(— the result 
follows. □ 

Lemma 3.3. Assume that Pi > ^2- For any y £ J there exist C3 < 00 
and 73 > such that 

sup Py {As ^ B for some s <T) < C3exp(— 73L) 

for L sufficiently large. 

Proof. We let sq = and, for A; > 1, define the following stopping times: 
tk = inf{t > Sk-i : At ^ {-L/2, 1/2^}, 
Sk = mf{t>tk:AteJ}, 
T = mf{t>0:At^B}. 
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Moreover, we denote by M{t) = supjA; >l-crk<t} the number of tips of 
arrows encountered by Ag by time t. Then for any site y £ J 

Py{As i B for some s < T) 

= Py{A„^ i B for some k < M{T)) 

< Py{A^, i B for some k < 2/3iT) + P{M{T) > 2piT) 

< Py{sk > T for some k < 2[5iT) + P{M{T) > 2(3iT) 

< 2PiTsupP,{si > r) + P{M{T) > 2piT). 

Since Ag has a drift toward J and the time between consecutive jumps has 
exponential bound, P{si > r) < C4exp(— 74L) for appropriate Ci< 00 and 
74 > (see the proof of Lemma 3.2). Furthermore, since EM(T) = PiT, 
large deviation estimates imply that there are C5 < 00 and 75 > such that 
P{M{T) > 2piT) < C5exp(-75r). □ 

Lemma 3.4. Assume that x E Bz, \\z\\ = 1 and Pi > (32- There exist Cq, 
76 G (0, 00) such that 

sup P^{At iB)< C6exp(-76L) 

for r and L sufficiently large. 

Proof. By decomposing according to whether Ag £ J for some s <T 
or not, we obtain 

Px{At <^B)< P^{As i J for ah s < T) + P{At ^B-A^eJ for some s < T). 

The first term on the right-hand side can be bounded using Lemma 3.2. For 
the second term, we first observe that 

P{At ^ B;As£ J for some s<T)< sup Py {As ^ B for some s<T) 

y&J 

and then apply Lemma 3.3. This completes the proof. □ 

Since there are {2L + 1)'^ sites in i?^, it follows from Lemma 3.4 and 
duality that there is a constant C7 < 00 independent of L such that for F 
and L sufficiently large 

P{r]Tix) = 2 for some x G B^) < J2 ^^(^T ^ B) 

< {2L + l^Ce exp(-76L) < CrL'^ < e/3. 



16 



N. LANCHIER AND C. NEUHAUSER 



Moreover, since each of the dual paths has a drift toward J, we can fix 
M > 0, say M = 3, so that for any e > 

P(any of the selected paths is not contained in 

\-ML, MLf at some time s<T)<s/3 

by choosing L sufficiently large. This shows that boxes that are sufficiently 
far apart are independent of each other with high probability and completes 
the proof of Proposition 3.1. 

4. Proof of Theorem 1. This section is devoted to the proof of part (b) of 
Theorem 1 which describes the behavior of the process for f3 > 1 and g close 
to in the generalist case a = (3. As already explained in the Introduction, 
if g = 0, then the process (t '■ — > {0, 1} reduces to the biased voter model 
with parameters 1 and /3 so that if /3 > 1, then P{(t{x) = 1) — > 1 as t — > oo. 
To prove that the pathogens still survive when 17 > is sufficiently small, 
we show that for M = 3 and any e > 0, we can choose L and F such that 
Proposition 3.1 holds for f3i = j3 and (32 = 1. Now that F and L are fixed 
and M = 3, we can assert that there exists (7c > small so that 

P(some infected host in [—ML,MLY' gives birth 

to a healthy host sent to B between time and time T) 

< {2ML + 1)^^(1 - exp(-5cr)) < e/3. 

This tells us that if (7 < (7c, then the set of occupied sites dominates the set 
of open sites in an oriented percolation process with parameter p=l — e. 
Here (z,n) occupied means that all sites in are occupied by pathogens 
at time nT. Finally, by taking e > sufficiently small so that percolation 
occurs with positive probability, Theorem 1 follows. 

5. Proof of Theorem 2. The proof of Theorem 2 also relies on a perturba- 
tion argument. In the case (3 = 00 and Ri < R2, the transition (1,0) (1, 1) 
is instantaneous, that is, unassociated hosts of type 1 become instanta- 
neously associated with a mutualist, provided that all the hosts of type 
1 are initially associated with a mutualist. The assumption Ri < R2 is to 
avoid the problem of births of isolated, unassociated hosts of the same type 
that are not accessible to mutualists. Under these assumptions, the process 

performs a biased voter model with parameters (3i = g and /?2 = 1- In 
particular, well-known results about the biased voter model imply that if 
g>l, then P{(lix) = 1) ^ 1 as t ^ 00. 

To extend the result to the region (3 > large, we prove that if hosts of 
type 1 become occupied by their associated mutualists quickly enough, then 
(,1 will evolve like a biased voter model in the space-time box B x [0,FL] 
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with probability close to 1. We first define on the same space as the 
biased voter model r]t introduced in Section 3 with Pi = g and /32 = 1 • At 
time T^'^, the host present at site w gives birth to an unassociated host of 
the same type which is then sent to x. At time U^''^ , the birth from w to x 
occurs only if the host at w is associated with a mutualist. To describe the 
evolution of the mutualists, we consider one more collection of independent 
Poisson processes, {Vn'^ ■n>l}, < ||x — w\\ < R2, with parameter /?. At 
time V^'"^ , we draw an arrow labeled with a 1 from it; to x to indicate that 
a mutualist (of type 1) present at site w gives birth to a mutualist at site x 
if this site is already occupied by a host of type 1 . We will prove that there 
exists Pj^"^ S (0,00) such that if /3 > PjJ^"^ and = 770 on B, then = 777- 
on Bz with ||2;|| = 1 at time T = TL with probability > 1 — e/3. Since boxes 
that are sufficiently far apart are independent of each other with probability 
close to 1, we can focus on [—ML,ML]'^ x [0,rL], M = 3, to estimate this 
event. Let x G [—ML, MVf and follow the line {x} x [0, TL] by going forward 
in time. Each time a host at w attempts to give birth at site x, we require 
that the next 1-arrow from w to x appears before the host at w is replaced 
or the host at x gives birth. A straightforward calculation shows that this 
event occurs with probability 

< min(rf ^ [/f for any y G Ml and 

If'"' < min(rf C/f for any y G A^) = + 2mY\ 

where m = gu^^ . Let us now denote by K{x, T) the number of unlabeled ar- 
rows and (5-arrows that point at site x by time T, and set Lm = [—ML, MVf. 
Then, by observing that ¥jK{x,T) = mT, and by decomposing the event to 
be estimated according to whether K[x,T) > 2mT or not, we finally obtain 

P{eT + VT on B,) < P{K{x,T) > 2mT) + 2mT ^ 

< {2MLy X {Csexp{--f8T) + 4m'^T{(3 + 2my^} 

for appropriate Cg < 00 and 73 > 0. Taking L and then f3 sufficiently large 
so that 

P{^^ ^rjT on B,)<e/3, 

and applying Proposition 3.1 imply that the set of occupied sites domi- 
nates the set of open sites in an oriented percolation process with parameter 
p = l — e. Here {z,n) occupied means that all sites in B^ are occupied by 
associated hosts of type 1 at time riT. This almost produces Theorem 2. 
Our last problem is that oriented site percolation has a positive density of 
unoccupied sites. To prove that there is an in-all-directions expanding region 
which is void of hosts of type 2, we apply a result from Durrett [7] which 
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Fig. 3. Process with nearest-neighbor interactions on the 200 x 200 torus at time 250. 
Left /3 = 2 and g = 0.5. The process starts with unassociated white hosts in J = (90, 110) 
and black hosts associated with a pathogen outside J . Right: /? = 4 and g = 2. The process 
starts with black hosts associated with a mutualist in J and unassociated white hosts outside 
J. In both pictures, gray sites refer to black hosts associated with a symbiont. 

shows that unoccupied sites do not percolate when e is close enough to 0. 
Since hosts of either type cannot appear spontaneously, once a region is void 
of one type, this type can only reappear in the region through invasion from 
the outside. This implies that P{^l{x) = 2) — > as t — > oo for any a; G Z"^ 
and completes the proof of Theorem 2. 

6. Proof of Theorem 3. This section is devoted to the proof of Theorem 
3. The method of the proof can also be applied to give a more explicit proof 
of Theorem 2 without too much more work. For the proof, we will assume 
that the birth rate A is not set equal to 1. In fact, we will prove Theorem 
3 first for A small and then change the time scale so that it holds for A = 1 
as well. We start by introducing the rightmost host-2 process r| and the 
leftmost symbiont-1 process £}, that is, 

rf = sup{x G Z : ^}{x) = 2} and = inf{x G Z : ^^{x) = 1}. 

We observe that because of nearest-neighbor interactions, — r| > 1 at any 
time t > 0. Moreover, 'd Gt = i] — — 1 denotes the number of sites between 
both processes, then changes as follows: 



for Gt = 0:r^ 




at rate A, 
at rate A^ 



and 



for Gt>0: 




at rate A 
at rate A 
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which suggests that drifts to the right if (7 < 1 (and to the left if (7 > 1). 
To make this argument precise, we wih prove the following lemma. 

Lemma 6.1. For g ^ I, there exists Eq > such that limr^oo 7^ x 
/o^l{G,=o} dt>eo. 

To deduce Theorem 3 from Lemma 6.1, we first observe that Er^ > A(l — 
g)£ot for 5 < 1- A large deviation estimate then implies that 

P(r2 > (1 - e)at) > 1 - Cge'^'^ 

where a = A(l — g)£o. This implies that if g <1, then — > 00 with probabil- 
ity 1. If g > 1, then Erf < —X{g — 1)^0^) and a similar argument implies that 
r| ^ — 00 with probability 1. The proof of Lemma 6.1 is based on a rescaling 
argument. The main objective is to prove that the fraction of time the host 
present at site rj + 1 is associated with a symbiont is greater than some 
positive constant. To be able to compare the particle system with oriented 
percolation process, we will artificially freeze by introducing the process 
seen from the interface, namely it{x) = + r| + 1). Prom this new point 
of view, the symbiont evolves on the half-line Z"*". To do this comparison, we 
need to extend some results on oriented percolation to oriented percolation 
in half-space, that is, the process in which all sites to the left of are closed. 

Oriented percolation in half-space. As explained above, our first objective 
is to investigate 1-dependent site percolation process in the half-space. Let 

= {[x ,n) : X -\- n is even and n > 0}. 

For any (x,n) G il, let Lv{x,n) define a random variable with the following 
property: If x < 0, then Lo{x,n) = 0, while if 2; > 0, then u;{x,n) = 1 with 
probability p and u;{x,n) = with probability 1 — p. The site {x,n) is said 
to be open (resp. closed) if a;(x,n) = 1 (resp. 0). Finally, 1-dependent with 
parameter p means that whenever {xi,ni), 1 < i < m, is a sequence with 
Xi>0 for any l<i<m where \\{xi,ni) — {xj,nj)\\oo > 1 for i / j, then 

P{uj{xi, Hi) = for 1 < i < m) = (1 - p)™. 

We say that {y,n) can be reached from (x,m) and write (x,m) — > {y,n) if 
there exists a sequence of points x = Xm, Xm+i, ■ ■ ■ ,Xn = y such that \xi+i — 
Xi\ = 1 for m<i<n—l, and uj{xi,i) = 1 for m < i < n. To formulate the 
next result, we let 

Wn = {y- ix,0) ^ {y,n) for some X € A} and = inf{n : = 0} 

for any A C 2Z+. Finally, if A = 2Z+, we will write Wn = and r = r^. 
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Lemma 6.2. If p> I — 3 then there exists ei > such that P{0 G 
W^2n) ^ £i for any n > 0. 

Proof. The proof closely follows ([5], Section 10) so we will try to 
be as brief as possible, and refer the reader to this reference for details. 
The first step is a slight modification of the contour argument applied 
to independent site percolation process. Let A = {0,2, ... , 2N} and C = 
{(?/,n): there exists xgA with (x,0) (y, n)}. We set 

D = {{a,b)eR^:\a\ + \b\<l} and W=[j{z + D). 

If the set C is finite, we denote by F the boundary of the unbounded compo- 
nent of (M X (—1, oo)) — W and orient the boundary in such a way that the 
segment from (2A^, —1) to {2N + 1,0) is oriented in the direction indicated. 
The boundary is a contour line, if it exists, starting at (2A'^, —1) and ending 
at (0,-1). There are at most 3™~^ contours of length m. Moreover, for a 
contour of length m to exist, there must be at least m/8 closed sites to the 
right of zero. To prove this point, we call a segment a line segment of the 
form X + F where x £ C and F is one of the sides of D. The site closest to 
the right of the segment is the site associated with the segment. Here, right 
and left are defined according to the orientation introduced above. We call 
segments of T which look like \, \ aiid y segments of types 1, 2, 3 and 
4, respectively. By construction, a site associated with a segment of type 1 or 
2 must be closed. Let rui and fhi be, respectively, the number of segments of 
type i and the number of segments of type i located on the left of zero. Since 
the contour starts at (2A^, —1) and ends at (0,-1), m-i +777-2 = 7773 + 7774 + 2 A^, 
so if the contour has length 777, then 7771 + r?72 >m/2. The same reasoning 
leads to mi + 7714 = 7772 + 7713 so that fh2 + fh^ <m/2. Now, since the sites 
located on the left of zero are closed with probability 1, we obtain 777-1 = 
and 777-2 = "T-s ^ m/4:, which implies (7771 + 7712) — {fhi + 777-2) ^ m/A. Finally, 
since a site in can be associated with at most two segments of type 1 
and 2, it follows that the number of sites on the right of zero that must 
be closed is > m/8. Noticing that the shortest possible contour has length 
2N + 4, one can conclude that 

00 

P(t[o.2^1 < 00) < J2 3'"(l-p)"^/8 = Cio(3(l-p)i/8)2^ 

m=2iV+4 

if p > 1 — 3~^ and the variables iv{x,n) are independent. To deal with the 
1-dependent case, we observe that there are nine sites in with (\m\ + 
|n|)/2 < 1, so for each F of length m there is a set of m/72 sites which are 
separated by more than one and which must be closed for the contour to 
exist. In conclusion, if p > 1 — 3"'^^, then P{t^ = 00) > 0. Finally, if we map 
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0, into itself by sending (x,m) i— > (x,n — m) and reverse the orientation of 
the graph, then 

P{W^ r\B^0)= P{W^ r\A^0). 

Taking A = 2Z+ and S = {0}, it then follows that P(0 G W2n) [ P{t° = 
oo) > 0. In conclusion, the lemma holds by setting ei = P{t^ = oo). □ 

The process seen from the interface. To prove Lemma 6.1 for g <1, we 
introduce the nearest-neighbor contact process in which a particle at x 
dies at rate 2A, gives birth at rate ^ur^ and sends its offspring to one of the 
neighbors at random in M^. The process is modified so that particles located 
in the interval {—oo,r^] are removed. More precisely, each time a particle in 

tries to give birth to a particle which is sent to a site in (— oo,r^], the 
birth is suppressed. Furthermore, if a particle is present at site r^ +1 when 
the process r| moves to the right, then this particle is removed. (The proof 
is similar in the case g > 1, provided one replaces ^~ by where particles 
give birth at rate fSvR^ but now die at rate 2\g.) The first step is to prove 
that dominates the process . 

Lemma 6.3. If g <1, the processes cmd can be defined on the same 
probability space in such a way that if the inclusion C holds at time 0, 
then C at any later time. 

Proof. Let us start by observing that if the rightmost process jumps 
to the right, a particle located at r| + 1, if it is present, is removed from both 
processes and and that, if it jumps to the left, both configurations stay 
unchanged. Therefore, it suffices to prove that the inclusion holds when r| 
is constant, say r| = at any time. This follows from a standard coupling 
argument so we just need to define a graphical representation that preserves 
the inclusion. We consider for any x, z>0 with |x — z| = 1 the independent 
Poisson processes {S^'^ : ?^ > 1} and {T^'^ : n > 1} with rate A^f and A(l — g), 
respectively. For any x, z>0, with < ||z — < R2, we also introduce the 
Poisson process {U^'^ :n > 1} with rate /?. The evolution of is as follows: 
At time S*^'^, the host present at site x gives birth to a healthy host of the 
same type, which is then sent to z. At time T^'^, the birth occurs only if the 
host at X is healthy. At time U^'^ , a pathogen (of type 1) at site x infects a 
host of type 1 at z if it is present. Finally, the process evolves according 
to the following rules: At time U^'^, a particle at site x gives birth to a 
new particle which is then sent to z. If the site is empty, the birth occurs. 
Otherwise, it is suppressed. At times T^'^ and 5"^'^, a particle present at z 
is removed. Such a coupling leads to the desired result. □ 
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We now introduce the process seen from the interface: r]t{x) = + '"t "I" 
1). To describe this process, we define the translation operators ri and t_i 
by setting 



[Tir]]{x) = r]{x + 1) and [r_i?7](x) 



r/(x — 1), if X > 1, 
0, if a; < 0. 



Then i]t evolves according to the following rules: 

1. A particle at x gives birth at rate (3i'R2 to a new particle which is then 
sent to a neighbor within the neighborhood J\f^. If the target site is empty 
the birth occurs; otherwise it is suppressed. Moreover, no births are allowed 
to the left of 0. 

2. Each particle dies at rate 2A. 

3. Depending on whether ^j^(r| + 1) = or 1, Tjt, respectively, shifts as 
follows: 

(Tim, at rate A, , (Tim, at rate A, 

Vt^ { J. ^ \ and Tjt^i J. ^ \ 

[T-irjt, at rate A, [T-ir]t, at rate Af?. 

To compare i]t with an oriented percolation process in half-space, we let 
r and L be two large integers to be fixed later, and scale space by setting 

B = [-L,L], ^{z) = Lz and B^ = ^{z) + B. 

Let Jz = ^{z) + (— L/5,-L/5), and say that the site {z,n) G is occupied if 
there is at least one particle in the interval Jz at time riFL. Let us fix z > 
even, and start the process with one particle in J^. The first step is to 
investigate the process with no deaths inside and modified so that any 
particle outside B^ is killed. We denote by fjt this new process. 

Lemma 6.4. Let 62 = 6"''^ and [5 > \. Then T and L can he chosen so 
that 

PinVL n Jz+1 0) > 1 - 2^2/3. 



Proof. A standard coupling argument implies that fit has more parti- 
cles if we increase the range of the interactions. So, we just need to prove 
the result when the offspring is sent to one of the two nearest neighbors. 
The idea is to prove that the rightmost particle in 1)1 reaches the right edge 
of Bz by time FL, and then stays inside Jz+i H Bz until time VL. Let 

Tt = sup{x G Bz : fit{x) = 1} and r = inf{t > : rj = $(z -|- 1)}. 

Then, on the set {r > t}, we have E(rj — tq) > (/? — X)t. In particular, since 
P > X, the parameters T and L can be chosen such that P{t > FL) < £2/3. 
This implies that the rightmost particle will reach Jz+i by time FL with 
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high probabihty. To prove that the rightmost particle does not leave Jz+i 
until time FL, we observe that 

rt ^ rt — 1 at rate < A and 

>/3, ifri<$(z + l), 



Tt ^ Tt + 1 at rate >, „ . ^ , . , ^ 

\0, if n = ^>(z + l). 

Then well-known estimates about random walks imply that 

P{3te[T,rL]:rtiJ,+i)<rLl^-^ < £2/3 

for L sufficiently large. This completes the proof. □ 

We now fix F and L such that Lemma 6.4 holds. To extend the result to 
the process ?7t, we just need to choose A > sufficiently small so that the 
probability a death occurs in the space-time region Bz x [0, TL] is smaller 
than £2/3. In other respects, since the result holds for the process modified 
so that any particle outside is killed, it follows that: If A > is small, 
then the set of occupied sites dominates the set of wet sites in an oriented 
percolation process in the half-space with parameter 1 — 6~^^. Lemma 6.2 
then implies that 

P(there is at least one particle in at time InFL) > £1 > 

for any integer n > 0, provided that ijq contains infinitely many particles. 
Now, it is easy to see that there exists a constant £3 > independent of n 
such that: If there is at least one particle in the interval Jz at time 2nrL, 
then the probability that £,i~{r^ + 1) = rjt{0) = 1 for at least one unit of time 
between times InFL and 2(n + \)rL is greater than £3. This tells us that 

1 

for some appropriate constant £0 > 0. Since Lemma 6.3 implies that {?/t(0) = 
1} C {Gt = 0}, Lemma 6.1 follows from the previous inequality. A time 
change now allows us to set A = 1, which then completes the proof of The- 
orem 3. 



7. Proof of Theorem 4. This section is devoted to the proof of Theorem 
4 which addresses coexistence of the symbionts in the neutral case g = l- To 
remind the reader, we assume that the symbionts evolve as previously but 
the hosts perform a threshold voter model according to the following rate 
at x: 

f 1, if car d{zGZ'':0< ||x-z|| <Ri 
i^j at rate I and {z)=j}>e, 

1 0, otherwise. 
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It is easy to see that the critical value for the infection rate, (k), is 
strictly bounded away from 0. Namely, if k = 1, then the symbionts per- 
form a basic contact process with death rate 1, provided 6 < ur^, and birth 
rate f5i'R2- Furthermore, since the contact process is monotone, /3^^'^(k) > 
for K > 1, from which our claim follows. 
To prove Theorem 4, we will compare the particle system viewed on suit- 
able length and time scales with a 1-dependent oriented percolation process 
in two dimensions. The properties of the process in the absence of symbionts 
was described in [7]. To apply his results, we introduce, for any x G Z*^, the 
house 

= [xiL, {xi + 1)L) X • • • X [xdL, {xd + 

where L is an integer to be fixed later and Xi denotes the ith coordinate 
of the vector x. We fix a <1/k such that 9 < criyR^, and say that Hx is 
good if it contains at least aL*^ hosts of each type. For x G Z*^, we define 
II 2; II 2 = (la^iP + ■ • ■ + la^dP)^^^ and set B2{x ,r) = {y: \\y — x\\2 < r} . We say 
that i?2(0,r) is good if for any x G 62(0, r) the house Hx is good. For z 
even for even n or z odd for odd n, we will say that {z,n) is occupied if the 
following two conditions hold: 

1. For any x £ B2{zKei,K), the house Hx is good at time nFL. 

2. For any i = 1, 2, . . . , k, B2{zKLei, KL) has at least one symbiont of type 
i at time nFL. 

Here, ei denotes the first unit vector, and K and F are large integers that 
will be fixed later. Note that the set B2{zKei, K) is defined on the rescaled 
lattice, whereas B2{zKLei, KL) is defined on the original lattice. We will 
prove the following result. 

Proposition 7.1. Let e>0 and 9 < vrJk. There exists Per G (0, 00) 

such that if (3 > I3^^\ then K , L and F can he chosen in such a way that 
the set of occupied sites dominates the set of open sites in a 1-dependent 
oriented percolation process with parameter 1 — e. 

The first step in proving Proposition 7.1 is to summarize the results of 
Durrett ([7], Section 2), which describe the behavior of the process in the 
absence of symbionts. To formulate the result we are interested in, we set 
Ri = L[Mi + 1) where L and Mi are large integers. 

Lemma 7.2 (Durrett). Let e > and 9 < vr-^/k. There exist Rq, Mq o,nd 
F such that the following holds: If Mi > Mq and B2{0,RoMi) is good at time 
0, then, for L large, B2{0,RqMi) is good until time TL and B2{0,2RqMi) 
is good at time TL with probability at least 1 — e/3. 
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The sets described in Lemma 7.2 provide an environment favorable to 
the survival of symbionts. To explain this, we introduce, for any type i = 
1,2, ... ,K, the processes t]1 defined by r]l{x) = 1 if ^t{x) = i and 7?j (x) = 
otherwise. Since a = 0, it is easy to see that, for i = 1, 2, . . . , k, the processes 
r]\ do not interact. We fix a type zG {1,2,. ..,«;}, and focus on the process 
vjl. The evolution of rfj- is as follows: 

1. Each particle dies at rate at most k and gives birth at rate (ivR^. 

2. A particle born at site x is sent to a site z chosen at random from 

Nl 

3. If the target site z is occupied by an unassociated host of type i, 
then the birth occurs. Otherwise, it is suppressed. 

The proof of Theorem 4 relies, like Theorems 1 and 2, on a perturbation 
argument. More precisely, we first prove Proposition 7.1 in the extreme case 
/3 = oo, and then extend the result to the region /? > large. We denote by 
ffj. the process vjl modified so that no births are allowed outside B2{0, KL). 

Lemma 7.3. Assume that B2{0,K) is good until time FL and that at 
time there exists x £ i?2(0, KL) with rjQ{x) = 1. // -R2 ^ ^dL and (3 = 00, 
then 

{x^lf: i}{x) =i}r\ ^2(0, KL) = {x£Z'^: f]i{x) = 1} for all t < FL. 

Proof. This is elementary geometry. To begin with, we cover the set 
-62(0, KL) with a finite number of Euclidean balls Bj, j G /; each of them has 
radius r = VdL. Then, it is easy to see that, for any j £ /, Bj contains at least 
one house. In particular, as long as B2{0,K) is good, Bj contains at least 
one host of type i provided that Bj C -82(0, KL). At any time < t < FL, let 
us pick one at random and denote by Xj{t) its spatial location. Now, since 
we have 

mm{\\Xj{t) - Xk{t)\\ ■.k^j}<R2 for all j £L and for ah t < FL. 

This implies that, for any x, z £ B2{0,KL) occupied by a host of type i, 
there exists a chain of sites xq = x,xi, . . . ,Xn = z such that the following 
two conditions hold: 

1. For A; = 1, 2, . . . , n, — < i?2- 

2. For k = 0,1, . . . ,n, the site Xk is occupied by a host of type i. 

In particular, since fjl starts with at least one particle in B2{0,KL) and 
B2{0,KL) is finite, all the hosts of type i are instantaneously invaded by a 
symbiont at time 0. It is easy to prove by induction that this holds until 
time FL. If a host of type i gives birth to an unassociated host which is 
sent to a site x G B2{0,KL) at time t, pick Xj{t) such that \\x — Xj{t)\\ < 
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i?2- Since Xj{t) is occupied by a symbiont of type i, the host at x will be 
instantaneously invaded. □ 

To extend the result to /3 > large, it is convenient to construct the 
process from a graphical representation. For any type i S {1, 2, . . . , k} and 
X E Z"', let {T^'^ :n > 1} be independent Poisson processes with rate 1. At 
time T^'^ the state of x flips to (z,0) if the set has at least 6 hosts of 
type i. For x e Z*^, let {C/^ :n > 1} be independent Poisson processes with 
rate p. At time U^, we choose at random a site z from AAJ. If a host of a 
certain type is present at site x, and a symbiont of the same type is present 
at site z, then the host at site x becomes associated if it is not already. So 
that Lemma 7.3 holds for /5 < co large, we now require the following two 
good events, denoted by Gi and G2, respectively: First, we need a quick 
invasion of the ball -62(0, KL) by the symbionts. More precisely, Gi will be 
the event that if at time there exists x € B2{0,KL) with rjf){x) = 1, then 
for all z G B2{0, KL), with z occupied by host i, the host present at site z 
becomes associated before another host attempts to give birth in -62(0, KL). 
To estimate P{Gi), we observe that, for any z € B2{0,KL), the host at site 
z can be reached in at most 4KL/R2 steps by a symbiont, that is, if the 
host at z is of type i, then there is a chain of sites xo,xi, . . . ,Xn = z with 
n < AKL/R2, satisfying the conditions 1 and 2 above and such that xq is 
occupied by a symbiont of type i at time 0. We denote by ukl the number 
of sites in the ball -62(0, KL). Then since the transition i — > j occurs at rate 
at most 1 and there are k hosts and vkl sites in B2{0,KL), new hosts are 
born at rate at most kurl- We set n equal to the integer part of ^KL/R2. 
Then 

P{Gi) > 1 - UKin 



> 1 - UKL 



KVKL + fi/n 
AKL kukl 

R2 HVKL + R2fi/^KL' 



Now that -62(0, KL) has been invaded, we secondly require it to remain fully 
occupied until time FL. In other words, G2 will be the event that given 
that at time all hosts are associated, each time a host is born at some 
site X S B2{0, KL), it becomes associated before another host is born in 
the ball B2{0, KL); this occurs from time to time FL. Let denote the 
number of times a host is born in -62(0, KL) from time to time FL. Since 
EA < kvklFL, we find that 

P[N > 2kvklFL) < On exp(-7iirL) 

for appropriate Cu < 00 and 711 > 0. If only one host in i?2(0, KL) is unasso- 
ciated, it becomes associated at rate at least j3. Births of hosts in ^2(0, KL) 
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occur at rate at most hvkl- Let X be a random variable with exponen- 
tial distribution with parameter f3 and let y be a random variable with 
exponential distribution with parameter ki^kl- Then 

P(G^) < P{N > 2KUKLrL) + P{Gl;N < 2KVKLrL) 

< Cu exp(-7iirX) + 2KUKLrLP{Y < X) 

< Cu exp(-7iiri) + 2K.UKLrL '^^^^ . 

[5 + kvkl 

The proof of Proposition 7.1 is now straightforward. Let e > and assume 
that i?2(0, A') is good and that .62(0, KL) has at least one symbiont of each 
type at time 0. Fix Ri = L{Mi + 1) and K = RqMi, then apply Lemma 7.2 
and choose L sufficiently large so that i?2(0, K) is good from time to time 
FL and i?2(0, 2K) is good at time FL with probability at least 1 — e/3. Now, 
increase L and then choose /? sufficiently large so that both probabilities 
and P{G2\Gi) are greater than 1 — e/3. To see that this produces 
the desired result, we observe that if -62(0,-?^) is good from times to FL, 
then Lemma 7.3 implies that, on Gi H (^2, the balls B2{—KLei, KL) and 
B2{KLei, KL) contain at least one symbiont of each type. This completes 
the proof. 

To deduce the existence of a nontrivial stationary measure /i from Propo- 
sition 7.1, we start the process ^4 from a product measure in which each 
host is associated with a symbiont and has density I/k. Then, we take the 
Cesaro average of the distributions from time to time T and extract a 
convergent subsequence. By Proposition 1.8 of [15], the limit /i is known to 
be an invariant measure. To see that fi has the desired property, we observe 
that if L is large, then the law of large numbers implies that {z, 0), z even, is 
occupied with probability close to 1. Moreover, if e > is small, well-known 
percolation results imply that, at any level n, the density of occupied sites 
is positive, which implies that ^{^'^{x) = i) ^0 for any i G {1,2, . . . , k}. At 
this point, we have proved that there is a critical value /3^^^ G (0, 00) such 
that if /3 > Pj^^, then coexistence occurs. 

To see that Pj,^^ can be chosen so that if /3 < (3j^^ then coexistence 
does not occur, we rely on a standard coupling argument. If we think of the 
process as being generated by the Poisson processes introduced above, it is 
easy to see that if Pi < (52, then the processes with parameters Pi and P2 can 
be defined on the same space, starting from the same initial configuration, 
in such a way that the process with parameter Pi has fewer symbionts of 
type i for any i £ {1,2, . . . , k}. This completes the proof of Theorem 4. 
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